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SUMMARY 


This study presents variational theorems for a theory of 
small motions superimposed on large static deformations and 
governing equations for prestressed beams on the basis of 3-D 
theory of elastodynamics . First, the principle of virtual work 
Is modified through Friedrichs *s transformation so as to describe 
the initial stress problem of elastodynaraics . Next, the modified 
principle together with a chosen displacement field is used to 
derive a set of 1-D macroscopic governing equations of pre- 
stressed beams. The resulting equations describe all the types 
of superimposed motions in elastic beams, and they Include all 
the effects of transverse shear and normal strains, and the 
rotatory inertia. The Instability of the governing equations Is 
discussed briefly. 


INTRODUCTION 


Small motions superimposed upon large static deformations 
have been tackled by a variety of investigators. And differential 
as well as variational formulations have been derived for both the 
so-called Initial stress and initial strain problems (see, e.g., 
refs. 1-3 j and references cited there). A classical variational 
formulation for the Initial stress problem Is deduced from a 
general principle of physics and has certain advantages over a 
differential formulation (see, e.g., ref. 3 > where the principle 
of virtual work is taken as fundamental). This yields only the 
stress equations of motion and the natural boundary conditions. 

The remaining equations of the Initial stress problem should be 
introduced as constraints, The constraints, however, can be 
removed through Friedrichs ’s transformation. This has been 
Illustrated by de Veubeke (ref, 4) for classical elastodynamics. 

All the past efforts reveal how the static and dynamic 
behavior of structures may significantly change by the presence 
of Initial stress or Initial strain. Among those, we mention here 
references 5-8 and references 9-12 on Initially stressed shells 
and plates, respectively. On Initially stressed beams, the works 
of Brunelle (ref. 13) and Sun (ref. l4) are cited. Brunelle 
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derived the governing equations for a prestressed, transversely 
Isotropic beam via the direct Integration of 3-D field equations. 
Sun studied the equations for a Timoshenko beam having an Initial, 
In-plane compressive stress by the use of both Trefftz's and 
Blot's formulations. 

The purpose of this Investigation is twofold. The first aim 
is to modify the principle of virtual work, and then to obtain a 
generalized variational theorem which describes an arbitrary state 
of Initial stress. The procedure used In achieving this Is 
analogous to the one used In reference 4, The second aim Is to 
construct the governing equations of anisotropic beams under 
initial stress by the use of the generalized variational theorem 
together with an incremental displacement field chosen a priori. 
The displacement field allows to include all the effects of 
transverse shear and normal strains, and the rotatory Inertia for 
the prestressed thick beam In which they are significant. The 
resulting equations describe all the types of superimposed 
ejxtensional, f lexural^ and torsional motions of thick anisotropic, 
elastic beam of uniform cross section. The dynamic Instability 
of the prestressed beam Is also discussed. 


SYMBOLS 


In a Euclldaen 3-space, Cartesian tensors are used, and 
Einstein’s summation convention Is Implied for all repeated Latin 
(1,2,3) and' Greek (1,2) Indices, unless Indices are put within 
parentheses . 

L, A; C length and cross-sectd-onal area of beam; Jordan curve 


which bounds A 


V, S 
S’, S" 


^1* ^a’ ^3 


u^, u“ 


m,n 

1 


P 


'’i 


e . 


10 


’lO 


0 


entire volume of beam and Its total boundary surface 

complementary subsurfaces of S, where stresses and 
displacements are_, respectively ^ prescribed 

a system of right-handed Cartesian convected coordi- 
nates; lateral coordinates and beam axis 

components of displac%ment vector, displacement 
functions of order (m,n) 

mass . density 

components of unit outward vector normal to S and C 
components of strain and symmetric stress tensors 
prescribed steady temperature field 
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mn 

ai= u^, 

rpHl j n 

ij 

■pin , n n ni ) n 

Qin,n pin,n 
(■•), ( ) 




(°), (*) 


Isothermal elastic stiffnesses and strain-temperature 
constants 

moment of Inertia of order (m,n) 

components of acceleration and traction vectors 

stress resultants of order (m,n) 

body force and acceleration resultants of order (m,n) 

effective load and external force of order (m,n) 

partial differentiation with respect to time, t, and x^ 

field quantities belong to the reference state and 
prescribed quantities 


C functions with derivatives of order up to and Including 

“ m and n with respect to space coordinates x^ and tlme,t 

FUNDAMENTAL EQUATIONS 


Consider a simply connected elastic body V+S, with Its 
boundary S, In a 3-D Euclidean space E. The elastic body Is 
referred to a x. -fixed system of Cartesian convected coordinates 
In this space. 'When this body Is prestressed, we distinguish two 
states of the body: Its reference (or Initial) and spatial (or 
final) state. The reference state Is considered to be self- 
equlllbratlng following static loading In the natural (or un- 
disturbed) state of the body at time, t=to. We may summarize, for 
ease of quick reference, the fundamental equations (see, e.g., 
ref. 2) In the form 


0 

ij,i 

+ p°fj 

= 0 

In V 




(1) 

s ® 

1 IJ 

-tj* 

- 0 

on S' 

> "1 - “1* 

= 0 on 

S" 

(2) 

0 

Ij 

^IJkl 

^kl 




In V 

(3) 


for this state. Here, p° Is the known mass density of the body 
material, s?. the symmetric stress tensor, f? the body force 
vector per unit mass In V, u? the displacement vector, n. the 
unit outward vector normal to S, u?* and t?* the prescrlDed 
displacement and traction vectors in the cimplementary sub- 
surfaces S" and S’ of S, e?. the linear strain tensor, and G. .v, 

(C. ., T = C.., , = C, T . . L the'^ Isothermal elastic stiffnesses. ^ 
xjkl jikl klio 

Now, suppose that an Infinitesimal (or small) motion Is 
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superimposed upon the reference state. For this motion, we have 
the following fundamental equations: 


(s..+s?u. ).+p°(f. - 

" 3-0 ir J ,r" ,1 ^ J 

a. ) 

J 

= 0 in V 

(4) 

n. (s . . + s ? u . ) - t* = 0 

ij xr j ,r' j 

on 

S ' 

(5) 

u. - u* = 0 on S" 

11 



(6) 

e . . = 1/2 (u. . + u_. .) in 

J 5 J J 3 

V 


(7) 

s . . = C . . , , ( e, T - 0a, T ) in 
xj xjkl kl kl 

V 


(8) 

u, - V* = 0 & u. - w* = 

XX X X 

0 

in V(to) 

(9) 


in the spatial stafe. In these equations, s.., u., t. and so on 
Indicate small Incremental quantities superlm^osei upon those of 
the reference state (i.e., s?., u“, t?). And a.=ii. is the 
acceleration vector , v* and * a^*e tfte prescrlfeed^dlsplacement 
vectors. 0 is an Incremental^prescrlbed steady temperature field 
and a..= a., the strain-temperature coefficients at constant 
stres^. Alii, V(to) is used to designate V at t=to. 

Equations (1)— (9) describe completely the initial stress 
problem of Interest. 


VARIATIONAL THEOREMS 


To begin with, we express a principle of Virtual work as the 
assertation 


/y(s®j + p“(f« + f^)6u^dV - p0a^6u^dV 

+ Jg(t“ + t*)6u^dS (10) 


in the spatial state. Here, y. . denotes the Lagrangian strain 
tensor, and it is given by 


Y 


ij 


= e , 


Ij 


+ 1/2 (u. 


u , 




( 11 ) 


In equation (10), through the use of equation (11), we first carry 
out the indicated variations, apply Green -Gauss integral 
transformations and combine the resulting surface and volume 
Integrals. Next, we recall the usual arguments on incremental 
field quantities (see, e.g., ref. 2), take into account equations 
(1) and (2), and finally arrive at the variational equation of the 
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form: 


with 


6J = (SJ = 0 
aa 


(12 a) 





+ s: u. 
ir j ,r 


+ fy p®(f^ - a^)6u^dV 


'^'^22 - ■*■ ®lr’^i),r^^i “ tljSUjdS (12 b) 

The variations of displacements are arbitrary and Independent in 
this equation. Hence, equation (12) leads evidently to the stress 
equations of motion (4) in V and the natural boundary conditions 
(5) on S, as the appropriate Euler equations. 


Variational Theorem: Let V+S denote a regular, finite region 

Df space (see, e.g., ref. 15) in H, with Its boundary S, and define 
the functional J whose the first variation Is given by equation 
(12). Then, of all the admissible displacement states u.eG^p, if 
and only If, the one which satisfies the stress equations or motion 


(4) and the natural boundary conditions (5) as the appropriate 
Suler equations, renders 6J =0. 


This is a one-field variational theorem In which equations (6) 
-(9.) of the Initial stress problem remain to be satisfied as 
jonstraints. 


To include the rest of equations of the Initial stress problem 
'.n the variational formulation, we introduce dislocation potentials 
ind use Priedrichs’s transformation, and we closely follow de 
^eubeke (ref. 4). Thus, we obtain the following theorem. 


Generalized Variational Theorem: Let V+S denote a regular, 

'inlte region of space In H, with Its boundary S (S’qS"=0 and 
)’uS"=S), and define the functional I whose first variation Is 
;iven by 


rith 


61 = 61^^ + 


' >S' - t|]6UjdS 

+ J3r!(u^ - uJ)6tj^dS 

^^22 ^ -^V ^®ij " ^ijkl^®kl " ©a^3_) 35£j_jdV 


(13 a) 


(13 b) 
(13 c) 
(13 d) 


hen, of all the admissible states of u^j^eG^p, 




■t . .and 
00 ^ 
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s. .eC^ , If and only If , those which satisfy the stress equations 
of^ motion (4) In V, the natural boundary Gondltlons (5) and (6) 
for displacements and tractions on S" and S’, the strain- 
displacement relations (7) in V, and the constitutive equations (8 
in V, as the appropriate Euler equations, render 61=0. 

In the generalized variational equation (13), the Incremental 
field quantities (s,., u. , t., and e,.) are varied independently. 
And this is a four-rleld^varlationai'^ theorem. The admissible 
states are not required to meet any of the fundamental equations o 
the initial stress problem but the Initial conditions (9) only. 


BEAMS UNDER INITIAL STRESS 


Geometry and Kinematics 

A straight elastic beam is embedded in the space E, The beam 
is of uniform cross section. A, and, it occupies a regular, finite 
region of Space ¥ with its boundary S in E. The total surface S 
consists of two right and left faces, A and A^, and a cylindrical 
lateral surface S-. The beam is referred to tne x. -system of 
Cartesian convected coordinates located at the centroid of A, . The 
x_-axls is chosen to be the beam axis, and the x -axes indicate 
the principal axes of A which is bounded by a Jo^^dan curve C. The 
beam is under an initial stress field in the reference state. 


The Incremental displacements of the prestressed elastic beam 
are taken of the form: 


(14) 


Here, the u^ are functions of x„ and time, t, only. These 
terms readily accommodate low-frequency extensional, flexural and 
torsional superimposed motions. However, it should be kept in 
mind that, in the case Of torsion, equation (l4) can represent 
only the displacements of beams of elliptic and circular cross- 
sections, and for all other sections, more terms should be retaine 
in the expansion. The displacement field (l4) is like the one 
Mlndlln (ref. l6) used in his recent derivation of the governing 
equations for a non-initlally stressed elastic bar. 


M=1 


u.(x.,t) = I 


m,n=0 


r m,n (m,n)-, 
[x^Xg uj_ ’ V3 


Stress and Load Resultants 
We define the stress resultants of order (m,n) : 


i(m,n) _ 


= /a 


(15 a: 
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This represents the weighted, averaged values of stress tensor 
over a cross section of the prestressed beam in the reference 
state . 

In addition, we Introduce the body force, aeGeleratlon and 
load resultants, and the moment of inertia of order (m,n): 


p(m,n) _ r m n f t*dA 

X A X ^ X X A X ^ X 


Imn = ! A ^>2 


(m,n) _ 


M=1 


= I 


^o(m,n) ^ ^ 


p,q=0 m+p,n+q 1 

r-p(m,n) po(m,n)-, ^ ^ xVv fs s<> 

^^1 ’ z -* 1 2^a > ai’ al-'^ 

M=1 


a(p,q) 


p,q=0 


•M=l 


[ (ppo(in+P-ljn+q) p qpo (m+p,n+q-l) ^^(p,q) 

p pQ (m+p,n+q)^(p,q) ^ 

3 i > 3 


joCm.n) j ^ + (np+mq)T? 

p,q=o 


+ + pT 

. „mo(m+p,n+q-l)-, (p,q) mo(m+p,n+q) (p,q) 

+ qT32^3 Ju. + T^^ u^^33 

. mo(ni+p,n+q-l)-, (p,q) •. 

^33,3 -'^1,3 ^ 

[ (pT«ta+P-l>n-^4) + qT.(m+P.n+q-l))^,jp.q) 

mo(rn+p,n+q) (p,q) 

^33 ^i,3 ■' 


0 (m+p-l,n+q) 
31,3 


o(m,n) _ 

p,q=o 


(15 b) 


Prestressed Beam Equations — Instability 

Now, we shall derive the prestressed beam equations by the 
se of the generalized variational theorem (13) together with the 
ncremental displacement field (l4). First, upon substituting 
he expansion (14) into equation (13 a), we find the variational 
quation (l6). In this equation, the variations 6u^*^ are 
rbitrary and independent for any choice of m(=0,l)^and n(=0,D, 
nd hence it evidently leads to the macroscopic equations of 
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motion (17) as follows: 


dx- = 0, m,n=0,l 

0 m,n=0 ^ ^ ^ 


( 16 ) 


n(m,n) ^ m(m,n) _(m-l,n) „(m,n-l) p(m,n) <m»n 

^1 31,3 “ 11 " '21 1 ^1 


+ p0(p{®^^) - = 0, m,n=0,l 

Here, Is the effective Initial load given by 


(17) 


Qo(m,n) ^ j^o(ra,n) ^ j^o(in,n) 


( 18 ) 


Similarly, we evaluate the variational equation (13 b) and 
obtain the natural displacement and traction boundary conditions 
In the form 


(m,n) _ ^*(m,n) 


= 0, m,n=0,l on S. 


^j(ra,n) ^ = 0, m,n=0,l on and A 

Here, S*=A^uAj^ and S"=S^, and n^^+l for A^ and n^^-l for A^. 

Upon using of equations (13 o) and (13 d) together with (14). 
we have the strain distribution: 


with 


= 11 ' i ='1^2 

m,n= 0 ,l ^ ' 


(m,n) ^ 

1 J 1 3 J J 3 1 


(20 a; 


+ (n+l)( 62 ^u(”-"-^l) + ] (20 b 


and the macroscopic constitutive equations: 


f \ M=1 

r j 

*^•1 -1 I 


T I ( £ ^ ^ ^ - a 0 ^ ^ ) (21 

IJkl I V+p,n+q^^kl “kl^ 

P > 4 


where we take the temperature Increment of the form: 
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( 22 ) 


M=i 

S(x.) = I 

^ m,n=0 ^ ^ ^ 

Lastly, the Initial conditions, based on equations (9) and 

( 14 ), 

u<»>n) . ^,(ni.n) , ^ -(m.n) _ ^ (23, 

jomplete the beam equations (cf., ref. 17, where non-lnltlally 
jtressed beams are treated) under an arbitrary state of Initial 
5 tress field. 

The beam equations of equilibrium may be derived similarly 
)n the basis of equations (l)-(3); they are not written out here ' 
.n order to conserve space. 

To examine the stability of the prestressed beam equations, 
re first consider the beam with a set of Initial forces x* Next, 
re replace x by a prescribed set x*' And, as usual, we arrive at 
L system of linear homogeneous differential equations which 
lescrlbes the Instability problem under consideration. The sets 
'.re defined by 

X = In L, In L, on A) 


X 




In L, In L, 


here X Is a monotonlcally Increasing factor, and whenever It 
eaches certain values the equilibrating reference configuration 
ecomes unstable. The behavior of the eigenvalues of this factor 
s to be Investigated In each particular case of Interest. Some 
xamples of Instability will be reported elsewhere. 
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